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Abstract 

We study the motion of pedestrians through an obscure tunnel where the lack of 
visibility hides the exits. Using a lattice model, we explore the effects of communica¬ 
tion on the effective transport properties of the crowd of pedestrians. More precisely, 
we study the effect of two thresholds on the structure of the effective nonlinear dif¬ 
fusion coefficient. One threshold models pedestrians’s communication efficiency in 
the dark, while the other one describes the tunnel capacity. Essentially, we note 
that if the evacuees show a maximum trust (leading to a fast communication), they 
tend to quickly find the exit and hence the collective action tends to prevent the 
occurrence of disasters. 

Resume 

Nous etudions la dynamique des movements de foules dans un tunnel dont la visibi¬ 
lity est tres reduite. Tout en particulier, nous nous interessons a des tunnels dont les 
sorties ne sont pas visibles. A l’aide de notre rnodele - un automate cellullaire - nous 
exploitons les effets de la communication inter-personnelle parmis les pietons sur la 
structure de la nonlinearite du coefficient de diffusion. Nous modelisons l’efficacite 
de la communication inter-personnelle ainsi que la capacite des sorties a l’aide de 
deux barrieres. 
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1 Introduction 


This Note deals with the following evacuation scenario: A possibly large group 
of pedestrians needs to evacuate a long and obscure tunnel. The lack of visi¬ 
bility is due to either an electricity breakdown or due to a dense smoke. The 
basic modeling assumption is that the pedestrians are equally fit, do not know 
each other, and also, are unaware of the precise geometry of the tunnel. We 
wish to build a lattice model to explore the effects of communication on the 
effective transport properties of these pedestrians. As modeling tool, we use a 
particular type of particle system, known as zero range process (abbreviated 
here ZRP), whose dynamics is affected by two thresholds. One threshold - 
called activation threshold - models pedestrians’s communication efficiency in 
the dark, while the other one - the saturation threshold - describes the tunnel 
capacity. From the modeling point of view, the activation threshold is open 
to many interpretations. In this Note, we associate the size of this threshold 
not only with the geometric level of the possibility of communication, but also 
with the willingness and ability of the pedestrians to process the transmit¬ 
ted information to make a decision towards orientation to a potential exit or 
choice of speed in the dark. We refer to this as level of trust. Essentially, a 
small activation threshold implies in this context a high level of trust. 

In Fig. 1 we sketch the meaning of the two thresholds, the precise mathemat¬ 
ical definition is given in Section 2. Each solid circle represents a pedestrian, 
whereas the associated open (bigger) circle represents its communication do¬ 
main and level of trust. On left bottom, the number of pedestrians in the cell 
is so small that their typical distance is larger than the radius of the commu¬ 
nication domain. On left top we see that if the number of pedestrian is large 
enough information can propagate throughout the cell. In other words, we as¬ 
sume that the information can be efficiently transmitted among the different 
pedestrians as soon as any single communication domain (the open circles) 
intersects at least one other pedestrian. Essentially, we need here a minimal 
degree of packing of the open circles, which is guaranteed in our scenario by 
the activation threshold. 

On the right bottom part we indicate that provided the number of pedestrians 
in the cell is smaller than the number that can be accommodated at bound¬ 
aries, then it increases proportionally to the number of pedestrians in the cell. 
If the number of pedestrians in the cell is too high (see right top), than the 
number of them exiting the cell per unit time saturates. 

It is worth mentioning that thresholds-biased dynamics have been discussed 
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Figure 1. Sketch of pedestrians moving through a cell of the obscure tunnel driven 
by a two-threshold biased dynamics. 

also for other transportation scenarios; compare e.g. [1,2,3] (group formation 
and cooperation in the dark) and [4,5] (collective dynamics of molecular mo¬ 
tors). This Note focusses on communication efficiency and is organized as 
follows: our transport model is described in Section 2, while Section 3 con¬ 
tains the hydrodynamic limit of the model as well as numerical illustrations 
exploring the effects of communication on the effective transport properties of 
the crowd of pedestrians traveling the obscure tunnel. 


2 The model 


We consider a positive integer L and define a zero range process (ZRP) [6,7] 
on the finite torus (periodic boundary conditions) A := {1,..., L} C Z. Fix 
N E Z + and consider the finite state or configuration space hi := {0,..., N} A . 
Given u = (cui,... ,col) £ hi the integer uj x is called number of particle at site 
x E A in the state or configuration c o. We pick A, S E {1,..., N} with S > A, 
the activation and saturation thresholds, respectively. We define the intensity 
function 


g(k) = { 


o 

i 

k-A + 1 


if k = 0 
if 1 < k < A 
if A < k < S 


( 1 ) 


[ S - A + 1 if k> S 


for each k E Z + . The ZRP we consider here is the Markov process uy G Q, 
with t > 0, such that each site x E A is updated with intensity g(co x (t)) and, 
once such a site x is chosen, a particle jumps with probability 1/2 to the 
neighboring right site x + 1 or with probability 1/2 to the neighboring left site 
x •— 1. For more details on this modeling strategy, see [6]. 


The intensity function relates to the hopping rates r( x,x±1 \LJ x (t)) = g(u x (t ))/2 
and coincides with the escape rate r( x,x ~ 1 'l(u x (t)) + A x,x+1 \u x (t)) = g(u x (t)) 
at which a particle leaves the site x. The thresholds intend to control the 
escape rate. Essentially, the activation threshold A keeps the escape rate low 
for all sites for which oj x {t) < A, regardless the number of particles on x. The 
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saturation threshold S holds the escape rate fixed to a maximum value for all 
sites for which c o x (t) > S, regardless, again, the number of particles on x. In 
the intermediate case, A < uj x (t) < S the escape rate increases proportionally 
to the actual number of particles on x, see (1). In the limiting case A = 1 and 
S — oo, the intensity function becomes g(k) = k, for k > 0 , and thus the well- 
known independent particle model is recovered. A different limiting situation 
appears when the intensity function is 1 for any k > 1 and 0 for k — 0. For this 
we find a ZRP whose configurations can be mapped to a simple exclusion- like 
model states (cf. e.g. [7]). Interestingly, we can tune between the two very 
different dynamics either by keeping S' = oo and varying A or by keeping 
A = 1 and varying S. 


3 Hydrodynamic limit. Threshold effects 


We study the hydrodynamic limit A”, L —> oo. Particularly, we exploit the fact 
that the intensity function is not decreasing and use well-established theories 
to derive the limiting non-linear diffusion coefficient and the limiting current in 
presence of the two thresholds. The Gibbs measure with fugacity z G M+ of the 
ZRP introduced above is the product measure v z (r}i)v z (r} 2 ) • • ■ v z (vl) on N A for 
any 77 = ( 771 ,..., r) L ) G N A with u z (0) = C z and v z {k) = C z z k /[g( 1) • • • g(k)] for 
k > 1, where C z is a normalization factor depending on z, A, and S, namely, 
1/C Z = 1 + SfcLx z k /[g(l) ■ ■ ■ g(k)]. To compute the mean value of the intensity 
function g , we use is z ( 0 ) with 77 ( 0 ) = 0 to get 


As a function of the activity, the expectation does not depend on the particular 
choice of the intensity function. 

For what concerns the hydrodynamic limit, a special role will be played by 
the density 

OO 

p(z) = J2 k "z(k). (3) 

k=0 

Note that p(z) is an increasing function of the fugacity, indeed, dp(z)/dz = 
[u z (rjl) — [y z {;pi)) 2 ]/z > 0. Hence, it is possible to define z(p) as the inverse 
function of p(z). We observe that p is defined for any positive z if A is finite 
and S = 00 . 

The evolution of the distribution of the particles on the space A under the ZRP 
with thresholds A and S can be described in the diffusive hydrodynamic limit 
via the time evolution of the density function p(x, t ) with the space variable x 
varying in the interval [0,1] and for any time t > 0; compare [ 6 ]. Consequently, 


v z [g (w*)] = J2 d( k ) = C zZ + C ' 


Jt-i 


-g{k -!) 


= z . 


( 2 ) 
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the continuous space density p(x,t ) is the solution of the partial differential 
equation 

Wt p = -T x J{e) with J{e) = ~\ D{p) T x p (4) 

where the macroscopic flux J(g) incorporates the effective diffusion coefficient 
D given by 

d 

D (p ) = -Q- p u m bM] • ( 5 ) 

The coefficient D is computed in terms of the mean of the intensity function 
evaluated against the single site Gibbs measure with fugaeity corresponding 
to the local value of the density. Consequently, D depends of the value of the 
thresholds. 


Let us discuss the effects induced by the two thresholds A and S on the 
diffusion coefficient. We first recall some known results which are valid in the 
two limiting situations A = 1 and S = oo ( independent particle model) and 
A = S ( simple exclusion-like model). In the first case, one has C z = exp{— z}. 
Hence, by (3), it holds p(z) = z. Recalling (2) and the definition of z, we 
have u Z (p) [g(ffii)] = v p [g(u;i)] = p. By using (5), the diffusion coefficient reads 
D(p) = 1. On the other hand, in the latter case, one has g(k) = 1 for any 
k > 1 and g(0) = 0. Hence, C z — 1 — z, and it holds p(z) = z/(l — z). Here, one 
birds the law D(p) = 1/(1+ p) 2 , cf. [8]. Now, we illustrate the general strategy 
to compute the diffusion coefficient D for arbitrary values of the thresholds A 
and S. To do so, we first compute p(z). The precise expression of the diffusion 
coefficient can be then obtained using the general recipe in equation (5) and 
recalling (2). Indeed, 


D(P) 






-1 


z=z(p) 


( 6 ) 


The explicit expression of dp(z)/dz appearing in ( 6 ) is lengthy, hence we omit 
it. 


Fig. 2 shows the behavior of the diffusion coefficient as a function of the local 
density, and parametrized by the values of the thresholds. The upper left panel 
refers to the case A = 1 and for different values of S: the simple exclusion-like 
model is recovered for S' = 1, while the independent particle model appears 
for S = oo. Similarly, the upper right panel illustrates the case with S — oo 
and for different values of A: here the independent particle model corresponds 
to A — 1 and the simple exclusion-like model is recovered for A = oo. In both 
the upper panels of the figure we see that for the independent particle case, 
the diffusion coefficient is constant with respect to the local density. 

Furthermore, in Fig. 2, we remark the loss of monotonicity of the function D(p) 
for values of p exceeding some critical value (depending on the thresholds A 
and S). The behavior of D(p) displayed in the lower left panel refers to the 
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Figure 2. Left panel, top row: Diffusion coefficient D(p ) vs. p for A = 1 and for 
different values of the saturation threshold, i.e., 5=1,2, 5, oo. Right panel, top row: 
Diffusion coefficient D(p) vs. p for 5 = oo and for different values of the activation 
threshold, i.e., A = 1,2, 5, oo. Left panel, bottom row: Diffusion coefficient D(p ) vs. 
p for A = 3 and for 5 = 3,4,10, oo. Right panel, bottom row: Diffusion coefficient 
D(p) vs. p for 5 = 10 and for A = 1 , 2, 5,10. 

case A = 3. Considering particularly the green curve corresponding to S — 10, 
one sees the onset of a double loss of monotonicity of the function D(p ): for 
small values of the density, D stays close to the simple exclusion-like behavior 
and decreases with p. After one first critical value of the density, it starts rising 
up, until it drops down again, when p exceeds an upper critical value. This 
indicates the presence of a double threshold for the intensity function given 

by (1). 

In the lower right panel of Fig. 2, the diffusion coefficient is plotted for dif¬ 
ferent values of A for the case S = 10. The red and the blue curves refer to 
the extreme independent-particle and simple-exclusion-like cases. When the 
activation threshold is varied two effects are prominent: the non-monotonicity 
of the diffusion coefficient with respect to density shows up and, at fixed local 
density, the diffusion coefficient decreases when A is increased. Interestingly, 
note that for local densities close to 8 this appears not to be true: increasing 
the activation threshold induces an increase in the diffusion coefficient. For 
this very particular regime the dynamics accelerates due to the increase in the 
activation threshold. In other words, it seems that close to the local density 
8 (which, depending on the population type, is close to the maximum pedes¬ 
trians density before asphyxiation starts off), higher mistrust speeds up the 
dynamics. 

Excepting the just mentioned non-intuitive regime, the effect of the two 
thresholds on the diffusion coefficient can be summarized as follows: the smaller 
the activation threshold A and/or the higher the saturation threshold S, the 
higher is the diffusion coefficient, and therefore, the quicker the dynamics. 
From the evacuation viewpoint, “a small activation threshold increases the 
diffusion coefficient” means that 
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higher trust among pedestrians improves communication in the dark 

and therefore the exits can be found easier. The size of the saturation threshold 
simply decides on the exit capacity. Consequently, a higher saturation thresh¬ 
old leads to an improved capacity of the exists (e.g. larger doors, or more exits 
[9]) and therefore the evacuation rate is correspondingly higher. 
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